Flow of hypermagnetic helicity in the embryo of a new phase in the electroweak phase 

transition 
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The dynamics of the magnetic helicity during the electroweak phase transition in the early Uni- 
verse is studied. It is shown that the boundary surface between symmetric (hypermagnetic) phase 
and Maxwellian phase with a broken symmetry is a membrana for the separation of the magnetic 
helicity. Assuming the total linking number of knots of hypermagnetic field is negative, it is proved 
that the helicity rising in the Maxwellian phase is left-handed. 
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The primordial magnetic fields in the early Universe 
could be sufficiently strong not influencing neither ex- 
pansion of universe nor Big Bang nucleosynthesis. If such 
fields survive under certain conditions after recombina- 
tion (z < 1100) they could be as seed fields in the galactic 
dynamo mechanism [l|, Q . 

An essential topological characteristic of the magnetic 
field which is global invariant in expanding universe with 
the changing structure of magnetic field at different scales 
is the magnetic helicity, for example, H = J cfixCB-A) in 
the case of the Maxwellian field. By a modern knowledge 
it can significantly affect the evolution of magnetic 
fields in galaxies. 

The cosmological magnetic field and its helicity can be 
formed as a result of phase transitions in the early Uni- 
verse and, in particular, in the electroweak phase tran- 
sition. In this phase transition the hypermagnetic field 
converts into the Maxwellian electromagnetic field. 

In this paper we try to study how the helicity of hy- 
permagnetic fields is related to the magnetic helicity of 
Maxwellian fields during electroweak phase transition. 
We show that during this phase transition on the sur- 
face separating the phases a separation of magnetic he- 
licity exists. The magnetic helicity being collected in 
Maxwellian phase in the course of this separation is con- 
served then in further expansion of the Universe and the 
subsequent formation of galaxies. As shown in 0,0; the 
phenomenon is due to neutrino asymmetry and parity 
violation (P-non-invariance) in weak interactions. In the 
absence of neutrino asymmetry in an ideal plasma the 
helicity is conserved. 

Let us consider a bubble (an embryo of the Maxwell 
phase) of the radius R, inside of the hot plasma in the 



early Universe at the moment of the electroweak phase 
transition with the temperature Tew ~ 100 GeV. Let 
us assume that this bubble is growing with the constant 
velocity, R{t) = v{t — tEw), where the velocity v itself 
(w = 0.1 — 1 according to [6|) is unessential (cancelled) in 
the solution to our problem. 

It is important for our calculation that the value 
{t — tEw) /tsw <C 1 is small, or that the tempera- 
ture during the phase transition remains constant at 
the moment Iew = Mq/2TI^ = 0.23 x IQ-i" c , 
where Mq = Mpi/l.QQ^/g* is given by the Plank mass 
Mpi = 1.2 X 10^^ GeV and by the degree of freedom 
g* ~ 100. This implies that the radius of the bubble is 
much less than the horizon size {2tEw — Ih = 1-44 cm), 
R-^Ih- More precisely, we shall assume that the radius 
of the bubble is much less than the scale of the mean 
hypermagnetic field, R <C r[Y /uy <C Ih- 

In the rest frame of the medium as a whole the in- 
duction equation governing hypermagnetic fields By = 
V X Y outside of the bubble reads: 



dBy 

dt 



= V X ayBy + 77yV^By, 



(1) 
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while the evolution equation for Maxwellian field B = 
V X A within the bubble is analogous, but with different 
value of the parameter a. This parameter characterizes 
P-non-invariance of weak interactions and it is a scalar, 
but not a pseudoscalar like the corresponding parameter 
< V • (V X v) > in the standard magnetohydrodynamics 
i. 

After the phase transition such a parameter of the he- 
licity is the following 

a = 2 X IQ-^GpTX-^ ^ c^^^U , (2) 

a 

while before the phase transition it takes the form 

ay = 3 X IQ-^g^a-^T ^ , (3) 
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where Gp '^s the Fermi constant, T is the temperature, A 
is the spatial parameter of the inhomogeneity of the neu- 

(A) 

trino gas; Ca = tO-5 is the axial constant of weak inter- 
actions (the upper sign for electron neutrino); g is the 
coupling constant for hypercharge field in the Weinberg- 
Salam model; ^j^^ = /i^^ /T is the dimensionless chemical 
potential of neutrino, a = e,/x,T, and coefficients of hy- 
permagnetic diffusion rj = jjy = (47r(7)~^ are given by the 
conductivity of plasma a 100 T and, practically, coin- 
cide for both phases. All coefficients ^i^^, ay, r]Y depend 
on the temperature (or time) by the Friedman law, but 
for given problem at the fixed temperature of the phase 
transition Tew: these parameters remain unaltered. 

Multiplying Eq. (m and its analogue for Maxwellian 
field by the corresponding vector potential and adding 
the analogous construction produced by evolution equa- 
tion governing the vector potential (multiplied by hyper- 
magnetic or magnetic field) after the integration over the 
space we get the evolution equation for the total helicity 
H = /(B • A)d^x + J (By • Y)d^x, where the integration 
is carried over the domains with Maxwell phase and hy- 
permagnetic field correspondingly. This equation takes 
the form: 



dH 
'dt 



-2 J {E-B)d^x- 

-^((E X A + AoB) •n)d25 + ..., (4) 



where the dots mean analogous terms for hypermagnetic 
field. We take into account the surface integrals which 
are omitted in problems for a monophase medium jT] as 
integrals over an infinite boundary of the domain. In 
our problem namely these integrals determine the flow of 
the helicity through the boundary of a bubble of the ra- 
dius R, on which a separation of the helicity takes place. 
Accounting for the boundary condition Af^ = cos9\yY^, 
where sin^ 9\y ~ 0.23 is the parameter of the standard 
Weinberg- Salam model, the integrals above are calcu- 
lated over the surface as following: 



dHy 
"dT 



— sm 



^(Ey X Y + YoBy)nyd^S, 



(5) 



where the unit normal vector ny = —Br = (—1,0,0) is 
directed inwards the bubble with Maxwellian phase (the 
phase with broken symmetry). 

The flow of hypermagnetic helicity density, penetrated 
inside the bubble through the surface at the moment of 
electroweak phase transition, is the pseudovector given 
by the formula 



S — nyhy{t) — ny 



1 



tEW 



{E'^ + B'^)/8tt is the energy density of the field. Here 
iirdR^ (t) is the volume of a thin spherical layer with the 
thickness d of the domain wall separating the two phases. 
For fields at the scale R, Y ^ ByR, the value of the flow 
is inversely to the thickness d , hy ~ d"^. 

It is not difficult to prove that the surface integrals are 
equal to zero, i.e. there is no separation of the helicity if 
we substitute the flat hypermagnetic field, Yq = Y^ — 0, 
Yx — Y{t)smkoz, Yy — Y(t) cos koz, (see estimates of 
the baryon asymmetry in [8|). But in the case of 3D- 
field with nonzero helicity the considered integrals are 
nontrivial. Let us consider the following potential of the 
hypermagnetic field with the number of the linked loops 
equal to n : 



Ye{t,p,e) 



-Y{t) cos 9 



(p2 + 1)2 

Y{t) sin 



(p2 + 1)2 

—Y{t)n sin 



1 + B{p-1)^ +b{p-lf 

p + cip-iy + 



(p2 + 1)2 

+iC + c)ip-lf], (7) 



where p = r/R, dependence on time is Y{t) — 2B[){t)/'KR 
(compare [10'|), and the coefficients b,c,B,C are calcu- 
lated below. We correct here the misprint in our paper 
■ 9,] where in the same Eq. ([7]) we missed the amplitude 
Y{t) for the component Yg{t,p,9). 

The hypermagnetic field By = V x Y near the surface 
of phase separation, < p — 1 ^ 1 , has the components 
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At the surface of bubble p = 1 our potential ([7]) and 
corresponding hypermagnetic field ([8]) are like in paper 
[lo| . Obviously, V • By = 0. We used the Lorentz gauge 
dY^/dx^ = to calculate the temporal component of 
hypercharge field 



This flow is analogous to the vector flow of the energy 
of a flat electromagnetic wave S = W^n, where W = 



3 



A straightforward calculation of the surface term ([5]) 
gives the following equation: 



dHvit) _ 27rsin2 6'w'n 
~dt ~ 3 



Rit)Y{t) 



Y{t 



T^dt , (9) 



where we substituted in the expression Ey 
VFq the gradient VYo, 



^dY/dt- 



1 



4 sin 9ea 



Y{t)dt 



4cos6'(l - 5p^)er 



(p2 



1 



(10) 



and took into account that in the case of the axial- 
symmetric configuration ([T]) the vector is indepen- 
dent of the coordinate 4>. The values in Eq. ^ including 
Vlo X Y = er{WYo)gY^ are calculated at the surface of 
bubble p — I. 

Hence the problem is reduced to the calculation Y{t) 
from the Faraday equation ([T]) which for the considered 
potential ([7]) and hypermagnetic field ^ at the boundary 
p = I can be rewritten by components as 



dBl 
dt 



ryy(V2By) 



' sin do 



'sm0B] 



aYY{t) cos ( 



dBY ,„9„ ^ ay d , 
— f - w V^By , = -^ — {tBD 
at r or ^ 



aYY{t) sin 6* 
2mt) 



[B + !)={), 



at r 



Q;yF(t)n sin 6* 



(C-2) 



(11) 



Let us pay attention to the zeroth result for Bj at the 
boundary p = 1 (see also in Eq. ([5])) that forces us 
to choose B — —1 in Eq. (fTTjl) for that component. 
The Laplacian in the l.h.s. of Faraday equation for the 
same component also vanishes at the boundary p = 1, 
{V^'BY)e = (5 ic)Y{t)nsine/2W'{t) = 0, if we choose 
c = -5/3. 

Then accounting for the other Laplacian compo- 
nents on the same boundary surface, the radial one, 
(V^By)^ = Y{t)ncose{2 - C)/R^{t) ^ 0, and the ze- 
roth (under conditions h = —2, B = —1) azimuthal com- 
ponent, (V^By)^ = Y{t) sme{3b - 3B + 3)/2R^{t) = 0, 
one gets from the first equation (jlip for radial compo- 
nent, 



Y R_2aY t]y{4.-2C) 

Y ^ R^ ^2 ' 



(12) 



while from the third equation (jlip for the azimuthal com- 
ponent we find 



Y R _ ayn(2 - C) 
Y^R^ R 



(13) 



Thus, in addition to the parameters B = —I, c ~ —5/3, 
b = —2 combining eqs. ([T2|) and (|13p we find the last 
parameter C{t) in Eq. 



Cit) 



2{n - n-i)ayi?-i + Ar^yR' 
nayR'^ + 2?7yi?-2 



(14) 



Substituting the parameter ((TH), e.g., into Eq. one 
obtains the ordinary differential equation for the ampli- 
tude Y{t), 



2al 



Y{t) R{t) nayRit) + 2r/i 



(15) 



In the realistic situation of finite conductivity a scale of 
the mean hypermagnetic field A — KrjY / ctY, where k > 1, 
should be much bigger than the diameter of the bubble 
in the new phase, i.e. the following inequality has to be 
satisfied: aYR{t) <C Kf^y. If a more stronger condition 
ayi?(i) <C 2riY/'n < kt/y is fulfilled, then from for 
the function By(i) = Y{t)/R{t) using ^ we get 



By(i) = Sy(t£;i4')exp 



-By(<£;H')exp 



63 



— ) {t — tEw) 

rjY 



0.001 



(16) 



where ByitEw) is the hypermagnetic field amplitude 
on the scale of the bubble, ay = ay (Tew), Vy — 
"riviTEw) are the constant coefficients at the moment of 
the phase transition, the sum = J^a^'^ai'^Ew) is the 
net neutrino asymmetry (neutrino degeneracy parame- 
ter); {t — tEw) /tRW ^ 1 is a small parameter for self- 
consistency of our problem (see above). 

Substituting the amplitude of the hypercharge field 
Y{t) = BY{t)R{t) on the surface of the phase separa- 
tion into the expression of the surface integral 
after the integration over time and division by the vol- 
ume of the spherical layer with the thickness d we get 
from ([6]) the value of the flow of hypermagnetic helicity 
density through the surface of the bubble, 



-3, 



/iy(t) _ 5 X 10 
G^cm d{cm) 



BYjtEw) 

1 G 



tEW ^ ' 



tEW 



(17) 



An unknown neutrino asymmetry at the moment of the 
phase transition is estimated by (^^/O.OOl) ~ 0.12. This 
estimate corresponds to the restrictions in Eq. (24) in 
paper obtained from the condition that the hyper- 
magnetic field survives against ohmic diffusion for spatial 
scales ~ rjY /aY- 



4 



Let us note that in order to avoid the screening of 
the hyperelectric field Ey and the temporal compo- 
nent Yq over the surface of the bubble, the thickness d 
of the domain wall should be less than the Debye ra- 
dius, d < td — \/ 3Tew 1 47re^rie ^ IO/Tew, that al- 
lows to estimate the factor in the formula ([T7| as 
d~^(cm) > 10^^/2. This means that a moderate hy- 
permagnetic field ByitEw) provides a huge flow of the 
hchcity density (fTT)) . 

Indeed, substituting into ([T7]) the value of hypermag- 
netic field at the moment of phase transition ByitEW^ 
estimated in Q as Byitsw) ^ 5 x 10^'' G, one gets 
h/G'^cm > Q.2^xl{)^'^[{t-tEw) /tEwf ■ Such huge value 
estimated at the moment of the growth of a bubble of the 
new phase, e.g, for R{t)/lH < [{t — tEw)/tEw] ~ 10"^, 
accounting for the following conservation of the net global 
helicity summed over different protogalactic scales, oc- 
curs much bigger than the helicity density of galactic 
magnetic field hgai ^ 10^^ G^cm, (see also estimates of 



the primordial magnetic helicity in paper 11|). 

We have to note that a growth of hypermagnetic 
field before the electroweak phase transition depends 
essentially (exponentially) on the neutrino asymmetry 
(Bvit) = exp[J^^{al.{t)/Ar]Y{t))dt] in a^-dynamo 
[5|). But in the expression for helicity (|17p a hyper- 
magnetic field -By(i) ~ ByitEw) is fixed at the moment 
of the phase transition, moreover, for a small bubble the 
answer is practically independent of neutrino asymmetry 
in the time-depended field i3y(i) given by ([16]). 

The single bubble of the Maxwellian phase inside of 
ambient symmetric phase with the potential given by 
Eq. ^ near the boundary, is a reasonable approxima- 
tion during the beginning of the phase transition before 
percolation (junction of bubbles). One can consider also 



another final step of the phase transition, when a new 
phase with broken symmetry prevails and a single bubble 
of the symmetric phase with hypermagnetic field inside 
exists. It is not hard to check that in this case the change 
of sign p— l>Otop— 1 <Oin the potential ([7]) gives 
the same components of hypermagnetic field inside the 
bubble p < 1. Let us note that in the considered ap- 
proximation ([7]) magnetic charges near the surface of the 
phase transition and over this surface itself are absent, 
V • By = 0. 

A choice of the negative sign of the helicity density 
p?)l if n < corresponds to the result 12| for the left- 
handed magnetic helicity in the same electroweak phase 
transition. That result is obtained for the mechanism of 
decay of linked loops of Z-strings leading to creation of 
magnetic monopole-antimonopole pairs at the ends of a 
decaying string, after which the reconnection of each such 
pair proceeds through junction by loops of Maxwellian 
field. 

Let us recall that the pseudoscalar n is the number of 
pairs of linking magnetic field loops entering the Gauss 
integral for magnetic helicity, H{t) = J d^xh{t, x) = 
n<I>i$2 0- This pseudoscalar changes the sign after one 
of the loops in a pair changes the sign (direction) of the 
flow 

For a single bubble of the symmetric phase the flow of 
the helicity density through the surface ([6]) preserves the 
value ([TT]) . Moreover, for the same n < this flow does 
not change the negative sign after the direction of the 
flow is changed, ny — ny = = (1,0,0). This well 
corresponds to the meaning of the problem: magnetic 
helicity of the Maxwellian field rises, unless helicity of 
the hypermagnetic field inside the bubble goes down. 
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